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Dijet cross sections are presented using photoproduction data obtained with the ZEUS 
detector during 1994. These measurements represent an extension of previous results, as 
the higher statistics allow cross sections to be measured at higher jet transverse energy 
(E 3 T et ). Jets are identified in the hadronic final state using three different algorithms, 
and the cross sections compared to complete next-to-leading order QCD calculations. 
Agreement with these calculations is seen for the pseudorapidity dependence of the direct 
photon events with E^ 1 > 6 GeV and of the resolved photon events with E 3 ^ 1 > 11 GeV. 
Calculated cross sections for resolved photon processes with 6 GeV < E J T et < 11 GeV lie 
below the data. 
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1 Introduction 



High energy collisions between photons and protons can produce jets in the final state. In 
leading order quantum chromodynamics (LO QCD), two types of processes lead to the pho- 
toproduction of jets. In direct processes the photon participates in the hard scatter via either 
boson-gluon fusion or QCD Compton scattering. In resolved processes the photon acts as a 
source of quarks and gluons, and only a fraction of its momentum participates in the hard 
scatter. This separation between direct and resolved photoproduction is only well defined in 
this way at leading order. To make a measurement which can be compared to calculations at 
any order, the variable x° BS is used to separate these two types of event |3J . The variable x° B 
is the fraction of the photon's momentum contributing to the production of the two highest 
transverse energy (Ej?*) jets. It is defined for the photoproduction of jets in positron-proton 
scattering as: 



7 " 2yE e 



where E e is the initial positron energy and rf et is the jet pseudorapidityQ. The inelasticity y is 

4(1- cos 9' e ) ' 



defined in the ZEUS frame as y = 1 — ^-(1 — cos 6' e ) where E' e and 6' e are the energy and polar 



angle of the outgoing positron. In a leading order calculation, direct processes have x^ = 1 



since all the photon's momentum participates in the production of the high transverse energy 



jets, while resolved processes have x° BS < 1 since part of the photon's momentum goes into 
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the photon remnant. Throughout the following, in both the data and the calculations, direct 
and resolved samples are defined in terms of a cut on x° BS rather than in terms of the LO 
diagrams. In a recent analysis by the HI collaboration, a similar variable was used to determine 
an effective parton distribution in the photon 0. 

In a previous analysis [|1[] dijet cross sections were measured using 1993 ZEUS data in the 
kinematic regime where the difference between the pseudorapidities of the two jets is small 
(|A?y| = \rf etl — rf et2 \ < 0.5). This condition constrains 8*, the angle between the jet-jet axis 
and the beam axis in the dijet centre of mass system, to be close to 90°. The cross section 
as a function of fj = (rf etl + rj : ' et2 )/2 then has maximal sensitivity to the parton distributions 
in both the photon and proton ||. In the comparison between data and Monte Carlo 
(MC) simulations based on the LO direct and resolved processes showed that the jet profiles, 
as described by the transverse energy flow around the jet axis, are poorly reproduced for jets 
with low E^ 1 produced in the forward (proton) direction. In the present analysis a comparison 
will be made with MC simulations which include multiparton interactions, and an improved 
description of the data is obtained. 

To compare data and theoretical cross sections based on next-to-leading order (NLO) QCD 
calculations, it is essential that similar jet definitions be employed for both the measurement 
and calculations. The dijet cross sections as a function of E^ and fj, for low and high x® BS , 
are measured in the hadronic final state using various jet definitions, including the algo- 
rithm. The resulting cross sections are compared to NLO QCD calculations at the parton level. 
The uncertainties due to hadronization effects are not yet theoretically estimated and are not 
considered in the comparison. After a brief description of the experimental setup, a discussion 
of the issues involved in the various jet definitions in both theory and experiment is presented, 
followed by our results and conclusions. 



x The pseudorapidity is defined as 77 = — ln(tan |) where is the polar angle with respect to the Z axis, which 
in the ZEUS coordinate system is defined to be the proton beam direction. 
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2 Experimental Setup 



In 1994 HERA provided 27.5 GeV positrons and 820 GeV protons colliding in 153 bunches. Ad- 
ditional unpaired positron and proton bunches circulated to allow monitoring of the background 
from beam-gas interactions. Events from empty beam crossings (that is bunches containing nei- 
ther positrons nor protons) were used to estimate the background from cosmic rays. The total 
integrated luminosity used in this analysis is 2.70 pb _1 with an uncertainty of ±1.5%. 

Details of the ZEUS detector have been described elsewhere The primary components 
used in this analysis are the central tracking system and the calorimeter. The central tracking 
system consists of a vertex detector || and a central tracking detector || enclosed in a 1.43 T 
solenoidal magnetic field. The uranium and scintillator calorimeter J7| covers 99.7% of the total 
solid angle and is subdivided into three parts: forward (FCAL) covering 4.3 > r] > 1.1, barrel 
(BCAL) covering the central region 1.1 > r\ > —0.75 and rear (RCAL) covering the backward 
region —0.75 > r] > —3.8, for a collision at the nominal interaction point. Each calorimeter 
part consists of an electromagnetic section followed by an hadronic section. The cells of these 
sections have inner face sizes of 5 x 20 cm 2 (10 x 20 cm 2 in the rear calorimeter) and 20 x 20 
cm 2 , respectively. A lead and scintillator calorimeter is used to measure the luminosity via the 
the detection of photons from the positron-proton bremsstrahlung process. This calorimeter 
is installed 107 m along the HERA tunnel from the interaction point in the positron direction 
and subtends a small angle at the interaction vertex ||. A fraction of the positrons scattered 
through small angles are detected in a similar lead and scintillator calorimeter positioned at 
Z = —35 m. 



3 Jet Algorithms 

Most of the previous measurements of jet cross sections at hadron-hadron colliders and in 
photoproduction at HERA have used some variation of a cone-based jet algorithm. In these 
algorithms, according to the standardisation of cone jet algorithms at the Snowmass meeting 
in 1990 0, jets consist of calorimeter cells (or, in a theoretical description, partons) i with a 
distance 

Ri = ~ V iet ) 2 + ~ ¥ et ) 2 < R (2) 

from the jet centre. Here fa and rji are the azimuthal angle and pseudorapidity of the cell (or 
parton), and R is the jet cone radius. In this analysis, the geometric centre of the cell is used 
to define the position. The parameters for the jet are calculated as: 

i 

V jet = -LE&W (3) 

i 

= iE^ 

E J T i 

in which the sums run over all calorimeter cells (or partons) belonging to the jet. Different 
approaches are possible to the choice of the 'seed' with which to begin jet finding, and to how 
and when overlapping jets are merged. The approach is not fixed by the Snowmass convention. 
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We use two different cone algorithms to determine dijet cross sections in photoproduction. The 
jet cone radius R = 1 is chosen for both algorithms. We also use a cluster algorithm, which 
does not suffer from these ambiguities. A further advantage of the cluster algorithm is that 
it is infrared safe to all orders, which is not always the case for cone algorithms [fLOR . In the 
following the three algorithms will be described in detail considering as an example the case of 
calorimeter cells. Identical algorithms are used in this analysis to define jets in the hadronic 
final state starting from the final state particles. 

In the first cone algorithm (EUCELL) a window in the r] — <fi space of the calorimeter cells 
is moved around to find those positions where the E? in the window is > 1 GeV to use as 
seeds. The jet quantities are initially calculated using the cells in a cone centred on the seed. 
Equations (0) and (|j) are then applied to choose the cells belonging to the jets and to update the 
jet quantities in an iterative procedure until a stable jet is found. Only the highest transverse 
energy jet is accepted, the cells within the jet are removed, and the whole process is repeated. 
In this way EUCELL produces no overlapping jets. 

The second cone algorithm (PUCELL) was adapted from the algorithm used by CDF jll]| and 



determines seeds by finding the single calorimeter cell of highest transverse energy and placing 
a cone around it. All the cells within the selected cone are assigned to this seed and excluded 
from the search for further seeds, which is then continued. The jet quantities are initially 
calculated using the cells in the seed and equations (|2|) and @ are then applied iteratively as 
for EUCELL until a stable jet is found. At this stage all jets are provisionally accepted. Thus 
it may happen that two stable jets overlap. If the overlapping transverse energy amounts to 
more than 75% of the smallest jet, they are merged, otherwise the overlapping energy is split 
such that cells are associated with the closest jet. 



In the cluster algorithm KTCLUS fl2 , [T3|| the quantity 



dij = ((Vi - Vj) 2 + ~ 0i) 2 ) min(£ ri , E T] f (4) 

is calculated for each pair of objects (where the initial objects are the calorimeter cells), and 
for each individual object: 

di = E 2 Ti . (5) 

If, of all the numbers [dij, di\, dk t i is the smallest then objects k and I are combined into a single 
new object. If however dk is the smallest, then object k is a jet and is removed from the sample. 
This is repeated until all objects are assigned to jets. As with the cone algorithms, eq. (|3]) is 
used to determine the parameters of the jets. It is also used to determine the parameters of 
the intermediate objects. 

Equations (0) and (|3]) imply that in a NLO calculation, two partons must be a distance 

r* = v^-^ + ^-^o 2 < J*:t E F s R ( 6 ) 

max(£/ Ti , rjTj) 

from each other to be combined, where E^ h is the transverse energy of parton k. This implies 
that if two partons have approximately equal transverse energy they may be separated from 
each other by as much as 2R and still satisfy eq.(^). However, as parton j does not then lie 
inside a cone of radius R around parton i and vice versa, one might with some justification 
also count the two partons separately. If one wishes to compare theory with measurement it is 
necessary to match the theoretical treatment of such cases to the operation of the jet finding and 
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jet merging criteria used experimentally. This is done by introducing an additional parameter, 
R sep , to the theory to restrict the maximum separation between two partons in a single jet 14 
Eq. (|) then becomes 



Rij < min 



-R. R 



m^(E Ti ,E Tj )^^ sep \- (?) 

The valid range of R sep is between 1R and 2R. For a NLO three parton final state, it is found 
2.0) • R corresponds to EUCELL, R sep = l-Rto PUCELL, and R = R, 



that R 



scp 



1.5 



'-scp 



to KTCLUS [15|. In this paper, all three jet definitions will be used for a comparison of the 
resulting dijet cross sections. An alternative approach would be to treat R sep as a parameter, 
and tune it in order to take into account possible theoretical uncertainties such as higher order 
contributions. However, in the present analysis this approach has not been followed and R sep 
is fixed by the functionality of each jet algorithm. 



4 Event Selection 

The ZEUS detector uses a three-level trigger system. The first level selects events used in 
this analysis with a coincidence of a regional or transverse energy sum in the calorimeter, 
and at least one track from the interaction point measured in the central tracking chamber. 
At the second level, at least 8 GeV total transverse energy, excluding the eight calorimeter 
towers immediately surrounding the forward beampipe, is required, and cuts on calorimeter 
energies and timing are used to suppress events caused by interactions between the proton 
beam and residual gas in the beam pipe [|16||. At the third level, a cone algorithm uses the 



calorimeter cell energies and positions to identify jets. Events are required to have at least 
two jets, each of which has E? LTJet > 3.5 GeV and v TLTJet < 2.0 or E? LTJet > 4.0 GeV and 
2.0 < rj TLTJet < 2.5. Additional tracking cuts were made to reject proton beam-gas interactions 
and cosmic ray events. 

Further cuts are applied offline. Charged current deep inelastic scattering events are re- 
jected by a cut on the missing transverse momentum measured in the calorimeter. To reject 
remaining beam-gas and cosmic ray background events, tighter cuts using the final Z— vertex 
position, other tracking information and timing information are applied. Two additional cuts 



are made [17 , based upon two different measurements of y: 



Events with a positron candidate in the uranium calorimeter are removed if y e < 0.7, 
where y e is the value of y as measured assuming the positron candidate is the scattered 
positron. 



A cut is made on the Jacquet-Blondel measurement of y fi~8f , yjB = J2i(Ei — E zi )/2E e , 
where E zi = Ei cos 9i, and Ei is the energy deposited in the calorimeter cell i which has 
a polar angle &i with respect to the measured Z- vertex of the event. The sum runs over 
all calorimeter cells. For any event where the scattered positron entered the uranium 
calorimeter and either is not identified or gives y e above 0.7, the value of yjs will be 
near to one. Proton beam-gas events will have low values of yjs- To further reduce the 
contamination from both these sources, it is required that 0.15 < yjs < 0.7. This range 
corresponds approximately to the true y range of 0.2 < y < 0.8. 
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These cuts restrict the range of the photon virtuality to less than ~ 4 GeV 2 , with a median of 
around 10~ 3 GeV 2 , which excludes deep inelastic scattering (DIS) events. 

To select dijet candidates with a particular jet algorithm, the algorithm is applied to the 
calorimeter cells. In each case, the jet transverse energy measured in the ZEUS detector is 
corrected as a function of rj CALJet and E^ ALJet . The variable E^ ALJet is used to denote the 
transverse energy of a jet before correction for the effects of inactive material. This correction 
is derived from the MC events described in the next section by comparing the true transverse 
energy of the jet, found by applying the algorithm to the final state particles, to the (lower) 
transverse energy measured in the calorimeter simulation, obtaining the average shift between 
the two transverse energies for each jet algorithm. The average shift in jet energies is around 
17% for all three jet algorithms, and varies between 10% and 25% depending upon rf Gt . The 
largest shifts occur at the boundaries between the FCAL and BCAL and between the BCAL 
and RCAL. No correction was applied to rj CALJet since, from MC, the average shift in rj between 
the particle and detector jets is less than ±0.05 for all i] values in the range used for the cross 
section measurements. 

The description of the calorimeter response to particles and jets in the MC has been tuned 



using several methods |19], p0| , including (i) the comparison of charged track momenta with 
calorimeter energy measurements, (ii) comparison of jet and positron variables in DIS events 
and (iii) the comparison of the measurement of the incident photon energy deduced from the 
energy of the positron measured in the small-angle positron calorimeter, to that calculated 
from energy deposits in the uranium calorimeter. The fivefold increase in statistics in 1994 
allowed the calorimeter energy scale to be studied in more detail than before, and these studies 
revealed a (6 ±3)% difference between data and MC. This difference was removed in the present 
analysis ( in the 1993 analysis [|l|, this difference was not corrected, but the possibility of such a 
discrepancy was allowed for in the systematic errors). Studies using jet photoproduction events 
allow us to assign an uncertainty of 5% to the calorimeter response for the jets used in this 
analysis P0|. Other studies have shown that of this 5% uncertainty, 3% arises from the absolute 



energy scale of the calorimeter flT9 |. 



After the jet energy correction, events are required to have at least two jets with > 6 GeV, 
— 1.375 < rf ei < 1.875, and \Ar]\ < 0.5. The MC gives a good description of the \Ar]\ distribu- 
tion around this region. For events with three or more jets, the two highest E^ 1 jets are used 
to calculate all jet-related event properties. This procedure is also employed later in all the 
theoretical and MC predictions shown. 

After these cuts about 25000 events remain, of which about 20% have x° BS > 0.75 (the 
exact number depending upon the jet algorithm). Events with a third jet which passes the Eij? 
and ?7 je * cuts comprise about 8% of the final sample. Of all events, 22% have their scattered 
positrons detected in the small-angle tagger, which is the fraction expected for a sample of 



photoproduction events |17| . No event from unpaired e or p bunches survives the selection cuts, 
implying that the non-ep background is negligible. The contamination from events with photon 
virtualities greater than 4 GeV 2 is estimated using simulated DIS events. This contribution is 
much smaller than the statistical errors and is therefore not subtracted from the data. 
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5 Comparison with Monte Carlo Simulation 



In Fig. ^ the x° BS distribution of the ZEUS data selected using the KTCLUS algorithm with 
(corrected) E J T et > 6 GeV, -1.375 < r) jet < 1.875, and \Ar]\ < 0.5 is shown (black dots). This 
x° BS is determined by using the corrected jet energies and corrected yjs- The correction to 
yjB is determined using MC generated events, by comparing yjB to the true y, as a function of 
the x° BS calculated using uncorrected variables. The peak at high x° BS due to direct photon 
processes and the rise at low x° BS due to resolved photon processes are both clearly visible. 
The sharp fall off for x° BS < 0.1 is a result of the E^ 1 and rf et kinematic cuts. 

The data are compared to the results of two LO QCD-based MC simulation programs, 
HERWIG58 ^ (solid line) and PYTHIA57 J22| (dashed line). All the MC events have been 
passed through a detailed simulation of the ZEUS detector and through the same jet energy 
correction procedure as was applied to the data. The GRV [^3] parton distributions are used 
for the photon and the MRSA |24| parton distributions are used for the proton. The simulation 
programs are based on LO QCD calculations for the hard scatter and include parton showering 
and hadronisation effects. The minimum transverse momentum of the partonic hard scatter 
(p^ in ) is set to 2.5 GeV in both HERWIG and PYTHIA. For both programs the direct and 
resolved photon processes are generated separately. 



as 



In the case of the resolved processes multiparton interactions (MI) are included [^, ^6 
an attempt to simulate the energy from additional softer scatters ('underlying event'), in both 
the dashed PYTHIA curve and the solid HERWIG curve. This has been shown to improve the 



simulation of the energy flow around the jet axis |27] 



In order to obtain the best agreement with the data the normalisations of the two processes 
were determined by allowing them to vary independently and fitting to the uncorrected x® B 
distribution. As a result, the cross section from HERWIG for resolved processes was scaled by 
1.2 with respect to the direct. The ratio of direct and resolved contributions using this scaled 
cross section was then 0.12, to be compared with 0.15 when using the unsealed cross sections 
within HERWIG. For PYTHIA the equivalent scale factor for the resolved cross section, and 
the ratio of direct and resolved, were found to be the same as for HERWIG within the precision 
quoted here. 

The dotted line shows the distribution for HERWIG without MI. For the MI, models based 
upon the independent statistical replication of scatters (eikonal models) are used which allow 
the generation of additional independent partonic scatters (with transverse momenta above 
pf n = 2.5 GeV for HERWIG and 1.4 GeV for PYTHIA) in resolved photon events. For HER- 
WIG the average number of scatters for events generated with these parameters is 1.05 and for 
PYTHIA it is 1.66. One effect of MI is to increase the number of events at low x^ B . However, 
even after the inclusion of MI with these parameters, the data still lie above the simulation 
at low x° BS . 

The uncorrected transverse energy flow around the jets is shown in Fig. |2|, for events in 
various bins of E^ ALJet and r] CALJet for KTCLUS jets, and is compared to the distributions 
from the HERWIG MC both with and without MI after full simulation of the detector. The 
jet profiles are described reasonably well by the MC with MI for most of the kinematic range, 
although there is still a tendency for MC jets to have too much energy inside the central region 
and too little energy outside this region, particularly for low E? jets in the forward region. This 
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Figure 1: x° BS distribution for KTCLUS jets with E J T et > 6 GeV, -1.375 < rf et < 1.875 and 
| At/ | < 0.5, where x° BS is calculated using corrected variables. The ZEUS 1994 data (black 
dots) are compared to the results of the HERWIG with MI (solid line) and without (dotted 
line) and PYTHIA with MI (dashed line) event generators after full detector simulation and 
scaling of MC cross sections (see text). The HERWIG cross section for resolved processes has 
been scaled by a factor of 1.2 with respect to the direct. The equivalent scaling for PYTHIA is 
the same. Only statistical errors are shown and in some cases are smaller than the black dots. 
The shaded area represents the direct process HERWIG MC events. 
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tendency is significantly stronger for MC samples which do not include multiparton interactions. 
However, we do not rule out the possibility that other models for the underlying event, or 
different MI parameters not investigated here, may provide a similar or better description of 
the data. 



6 Resolution and Systematics 

The resolution of the kinematic variables has been studied by comparing, in the MC simulation, 
jets reconstructed from final state particles (hadron jets) with jets reconstructed from the 
energies measured in the calorimeter (detector jets), and by comparing y JB with the true y. 

The distribution of the difference in fj between the hadron and detector jets has a mean of 
zero, a width of 0.15 units and depends weakly on fj, exhibiting shifts of about 0.05 units close 
to the boundaries between the BCAL and the FCAL or RCAL. The resolution in x° BS is 8% 
at x° BS = 0.75. For and y, the resolutions are 15% and 0.09 units, respectively. 

The jet cross sections presented in this analysis refer to jets in the hadronic final state. The 
MC samples have been used to correct the data for the inefficiencies of the trigger and selection 
cuts and for migrations caused by detector effects. The correction factors are calculated as 
the ratio N true /N rec in each bin. N tTue is the number of events generated in the bin and N rec 
is the number of events reconstructed in the bin after detector smearing and all experimental 
cuts. The final bin-by-bin correction factors are between 0.5 and 1.5 for all the cross sections 
measured. The dominant effect arises from migrations over the E^ 1 threshold. 

The sensitivity of the measured cross sections to the selection cuts has been investigated by 
varying the cuts on the reconstructed variables in the data and HERWIG MC samples and re- 



evaluating the cross sections f20|| . In addition, the cross sections were re-evaluated using a ratio 
of the direct and resolved contributions derived from the cross sections from HERWIG without 
additional scaling (direct/resolved=0.15), and by using the PYTHIA sample. They were also 
evaluated using the HERWIG model with and without multiparton interactions. These effects 
are included as systematic errors on the cross sections, and are correlated to some extent. The 
possibility that the detector simulation may incorrectly simulate the detector energy response 
by up to ±5% has also been considered, as mentioned in section (|. This effect is added 
in quadrature to the overall normalisation error of 1.5% arising from the uncertainty in the 
measurement of the integrated luminosity. This principal correlated uncertainty is indicated in 
the figures as a shaded band and should be added to the other systematic errors to give the 
overall uncertainty. 



7 Results 

The measured cross sections are now discussed and compared to theoretical expectations. The 
cross section is first measured over the whole x° BS region and its shape is compared to that 
of MC expectations. This cross section includes x° BS values down to 0.05, the lowest value 
allowed by the other kinematic cuts. At the lower values of x? f BS ' , the jet profiles and Fig. p] 
indicate discrepancies between the data and the MC simulations. Nevertheless, this cross 
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Figure 2: Uncorrected transverse energy flow (l/N)dE T /d8rj around the jet axis, for cells within 
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section remains interesting as its shape is less biased by kinematic cuts than those of the cross 
sections to be discussed in section 7.2. We compare the shape to MC simulations which include 
models for MI, parton showering and hadronisation, but have large scale dependences due to 
the fact that they include only LO matrix elements. 

Next, x® BS cuts are applied to select regions where contributions arising from an underlying 
event - which may be responsible for the low-x^ discrepancy in Fig. |T]- are reduced and hence 
NLO QCD can be expected to provide a better description of the jet production process. The 
cross sections measured here in the hadronic final state are compared to NLO QCD calculations 
of partonic cross sections. These calculations have a reduced scale dependence but do not 
include parton showering beyond a single branching. MI and hadronization effects are also not 
included since no theoretical estimation of these two contributions is yet available for these 
calculations. This uncertainty is not considered in the following comparisons. 



7.1 Cross Sections without x° BS Cuts 

The cross section da/dfj for ep — > e+ dijets +X in the range \Arj\ < 0.5, 0.2 < y < 0.8 and for 
virtualities of the exchanged photon less than 4 GeV 2 is shown in Fig. |3] and given in table |1] 
for the KTCLUS algorithm, requiring E° T et > 6 GeV. 

The cross section rises from around 0.2 nb per unit of pseudorapidity at fj = — 1 to around 
3 nb per unit of pseudorapidity for fj > 0.25. The data may be compared with the predictions 
of the HERWIG MC using the direct /resolved ratio of 0.15 given by HERWIG. While the 
simulation can describe the shape of the cross section, these predictions fail to describe the 
overall normalisation, requiring an ad hoc multiplicative scale factor of about 1.8 to agree with 
the data. Such a factor is not unreasonable given the scale dependence of the MC. Fig. |] shows 
various predictions of the HERWIG MC after including the factor of 1.8. With the value of 
pram _ 2 5 GeV used here, the data slightly favour the GRV parton distribution |23[] with MI. 



The LAC1 [ESI or the GRV distribution without MI also gives reasonable description of the 



data. However the LAC1 distribution with MI is ruled out. 

The effect of MI in the simulations is a strong function of the choice of the photon parton 
distributions, in particular the gluon component, which is where the major difference between 
LAC1 and GRV lies. Additionally, it should be noted that the effect of MI is also a strong 
function of the choice of p™ m f26f . No comparison is presented here with NLO perturbative 
QCD calculations since they do not include MI. These comparisons are performed in the next 
subsection, after applying x° BS cuts to reduce such effects. 

7.2 Cross Sections with x° BS Cuts 
Two regions have been selected: 

1. x° BS > 0.75: direct photoproduction. 

2. 0.3 < x° BS < 0.75: resolved photoproduction excluding the \ow-x^ BS region. 



10 



ZEUS 1994 



2 - 



mm 



T 



= 6 GeV, no cut on x_P B s 



y 



• ZEUS (KTCLUS) 

HERWIG (xl.8): 

GRV(LO) 

GRV(LO)[no MI] 
LAC1 

LAC 1 [no MI] 




principal correlated error 



-0.5 



0.5 



1.5 



*1 



Figure 3: da/dfj for ep — > e+ dijets +X in the range \Ar}\ < 0.5, 0.2 < y < 0.8 and for 
virtualities of the exchanged photon < 4 GeV 2 and for integrated above E™ m = 6 GeV. The 
cross section is measured using the KTCLUS algorithm and is compared to the expectations of 
various HERWIG MC simulations (see text). The errors bars represent the combined systematic 
and statistical uncertainty, excluding the principal correlated uncertainties which are shown in 
the shaded band, see text. 
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For each x° BS region, the cross sections da/dfj for ep — > e+ dijets +X in the range \Arj\ < 0.5, 
0.2 < y < 0.8 and for virtualities of the exchanged photon < 4 GeV 2 are measured for four 
different values of the E^ 1 threshold, E™ m = 6, 8, 11 and 15 GeV. The cross sections are 
measured for the three different jet algorithms discussed in section 3. The results are given in 
Tables ^| to [7| and are displayed in Fig. |] together with the results of the NLO QCD calculation 



from Klasen and Kramer [29] using the NLO GRV [23] parton distributions for the photon and 
the CTEQ3M |3l| parton distributions for the proton and employing two different values of 
the i? sep parameter: R sep = 1 (solid curve) and R sep = 2 (dashed curve). Since the jets may 
be accompanied by other soft gluons (outside the jets), there is a potential problem when the 
two jets have the same E^ 1 . The infrared singularity associated with summing the soft gluon 
contributions is usually cancelled by the singularity coming from the one-loop contribution. 
For two jets with the same E^, some of the phase space for the soft gluon terms is restricted 
and an incomplete cancellation may occur in some calculations. As a consequence, Klasen and 



Kramer |29| have allowed the second jet to have an E 3 ^ less than E™ m if the third (unobserved 
gluon) jet has a transverse energy of less than 1 GeV. However, the cross section is then sensitive 



to changes in the value used for this cut on the third jet. Harris and Owens |JT]] have applied 
a low cutoff on the energy of the very soft gluons and found that the dependence of the cross 
sections on the value of the low energy cutoff used is much less than the quoted errors on the 
data. These different approaches account for the differences between the theory curves shown 
later. 

The photoproduction cross section for x° BS > 0.75 and Ejt %n = 6 GeV (Fig. |]a) rises 
from around 0.2 nb at fj = — 1 to a maximum value of around 1.8 nb at fj = 0, decreasing 
back to 0.2 nb by fj = 1. This decrease arises from the cutoff on the minimum E^ 1 and 
the cuts on y. The EUCELL jet cross sections are systematically higher than the PUCELL 
cross sections, which in turn are slightly above the KTCLUS cross sections. This behaviour is 
qualitatively similar for the higher E™ m cross sections (Figs. f|b-d), where the maximum value 
of the cross section falls and occurs at steadily higher fj as the minimum Ej? cut increases. The 
PUCELL and KTCLUS cross sections are in good agreement with the NLO curve calculated 
with R sep = 1 for all fj and for all four E^ thresholds, except for the most negative values of 
fj in the lower E™ m cross sections, where the trend is for the calculation to lie above the data. 
The -R sep = 2 curve lies above all the data at most values of fj. In the data the separation 
between EUCELL, PUCELL and KTCLUS becomes less significant at higher Ej?. However, 
the separation between the two theory curves remains significant. 

The photoproduction cross section for 0.3 < x® BS < 0.75 and E™ in = 6 GeV (Fig. |]e) rises 
from around 0.8 nb at fj = —0.25 to a maximum value of 1.5 nb for PUCELL and KTCLUS, 
and of 3 nb for EUCELL, at fj — 0, followed by a decrease back to 0.2 nb by fj = 1.5. The 
EUCELL jet cross sections are again systematically higher than the PUCELL cross sections 
which are again slightly above those for KTCLUS. This behaviour is once more qualitatively 
similar for the higher E^ cross sections (Figs. f|f-h), where the maximum value of the cross 
section falls and occurs at steadily higher fj. In the data the separation between EUCELL and 
the two other jet algorithms is larger than in the direct factor of two at the lowest 

values - but again becomes less significant at higher E^. In the theory, the differences between 
the curves with different R sep again show the same trend as the data, with the R sep = 2 curves 
being higher than those for i? sep = 1. However, for the cross sections with E™ n = 6 GeV 
and 8 GeV, the NLO QCD curves lie significantly below the data. For higher E™ m values the 
calculations are broadly consistent with the data. 
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Figure 4: da/dfj for ep — > e+ dijets +X in the range \Arj\ < 0.5, 0.2 < y < 0.8 and for 
virtualities of the exchanged photon < 4 GeV 2 and for E™ m = 6, 8, 11 and 15 GeV. Figures 
(a-d) are the cross sections measured in the range x° BS > 0.75; figures (e-h) are for the range 
0.3 < x° BS < 0.75. The cross sections are measured using three different jet algorithms and 
are compared to NLO QCD calculations using R sep = 1 (solid curves) and R sep = 2 (dashed 
curves), see text for details. The errors bars represent the combined systematic and statistical 
uncertainty, excluding the principal correlated uncertainties, which are shown in the shaded 
band (see text). The band indicates the maximum uncertainty for the three jet finders. The 
individual uncertainty for each jet finder is given in the table. 
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In Fig. |5] the KTCLUS jet cross sections are shown again, with Klasen and Kramer's NLO 
QCD calculations (with R sep = 1) employing two different parton distribution functions for the 
photon - namely those of NLO GRV(solid curves p3| ) and GS (dashed curves P2| ). It can 
be seen that the agreement is in general good for both distribution functions, except in the 
two lowest Ej? regions of the resolved cross section (Fig |5|e,f ) where the QCD calculations are 
significantly below the data. Perhaps surprisingly, the difference between the photon parton 
distributions is largest in the direct photoproduction region. This is due to differences between 
the quark distributions in the photon for x 7 > 0.8, where they are poorly constrained by photon 
structure function measurements at e + e~ colliders. These differences persist at high E 3 ^ . Also 
shown in Fig. [|is a NLO QCD calculation (again with R sep = 1) from Harris and Owens using 
NLO GRV MS for the photon, and NLO CTEQ4 MS for the proton. At high x° BS there is 
again good agreement with the measurements, but at low x° B the disagreement in the two 
lowest Ej 1 regions is large. At higher Ej? values, the data and calculations are consistent. 



8 Conclusions 

Photoproduced dijet cross sections da/dfj have been measured in the hadronic final state for 
different kinematic regions and are found to be consistent with the general expectations of 
QCD, in the sense that both resolved and direct processes are observed in the data. 

Quantitatively, it is found that Monte Carlo simulations both with and without multiparton 
interactions are capable of describing the fj dependence of the cross section when no x° BS cuts 
are applied, although simulations which use multiparton interactions to simulate an underlying 
event are slightly favoured and also give a better description of the jet profiles. 

The measured cross sections vary by up to a factor of two when different cone or clustering 
algorithms are used for the definition of jets. This behaviour is similar to that predicted from 
the theoretical calculations by choosing the R sep parameter in order to reproduce the different 
jet algorithms. 

Comparison of the direct photon cross sections (x° BS > 0.75) with NLO QCD calculations 
shows good agreement in both shape and magnitude over a wide range of E^ 1 and rf et and 
for the three different jet definitions. It also displays a sensitivity to the photon structure at 
large x T 

Calculations for the resolved photon cross sections in the region 0.3 < x° BS < 0.75 which 
include jets with 6 GeV< E 3 T et < 11 GeV are found to lie below the data. However, for higher 
jet energies the calculations are consistent with the data. 
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Figure 5: da/dfj for — > e+ dijets +X in the range |At/| < 0.5, 0.2 < y < 0.8 and for 
virtualities of the exchanged photon < 4 GeV 2 and for E™ m = 6, 8, 11 and 15 GeV. Figures 
(a-d) are the cross sections measured in the range x° BS > 0.75; figures (e-h) are for the range 
0.3 < x? t BS < 0.75. The cross sections are measured using the KTCLUS jet algorithm and 
are compared to NLO QCD calculations from Klasen and Kramer, for two different parton 
distributions in the photon , GRV and GS, and from Harris and Owens using GRV (see text for 
details). All three calculations use R sep = 1. The errors bars represent the combined systematic 
and statistical uncertainty, excluding the principal correlated uncertainties which are shown in 
the shaded band (see text). 
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Table 1: The cross sections for KTCLUS for the whole x° range. The third and fourth 
columns represent the statistical and systematic uncertainties, excluding the principal corre- 
lated uncertainties which are shown in the fifth column, see text. 
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Table 2: The cross sections for PUCELL and x° BS > 0.75. The third and fourth columns 
represent the statistical and systematic uncertainties, excluding the principal correlated uncer- 
tainties which are shown in the fifth column, see text. 
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Table 3: The cross sections for PUCELL and 0.30 < x° BS < 0.75. The third and fourth 
columns represent the statistical and systematic uncertainties, excluding the principal corre- 
lated uncertainties which are shown in the fifth column, see text. 
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Table 4: The cross sections for EUCELL and x° BS > U.75. The third and fourth columns 
represent the statistical and systematic uncertainties, excluding the principal correlated uncer- 
tainties which are shown in the fifth column, see text. 
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Table 5: The cross sections for EUCELL and 0.30 < x° BS < 0.75. The third and fourth 
columns represent the statistical and systematic uncertainties, excluding the principal corre- 
lated uncertainties which are shown in the fifth column, see text. 
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Table 6: The cross sections for KTCLUS and x° BS > U.75. The third and fourth columns 
represent the statistical and systematic uncertainties, excluding the principal correlated uncer- 
tainties which are shown in the fifth column, see text. 
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Table 7: The cross sections for KTCLUS and 0.30 < x° BS < 0.75. The third and fourth 
columns represent the statistical and systematic uncertainties, excluding the principal corre- 
lated uncertainties which are shown in the fifth column, see text. 
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